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I. INTRODUCTON AND MOTIVATION 

Classical solutions of Yang-Mills theory can provide 
new insights into non-perturbative physics, giving hints 
on infrared properties (related to confinement or chiral 
symmetry breaking, for example). Among the known 
classical solutions, we can name several configurations 
which have been investigated and used in studies of 
hadronic phenomenology: the instanton, meron, dyon 
and monopole (for a review, see [l[). Furthermore, in 
spite of the dyon being a kind of generalized monopole so- 
lution |2| , and the existence of a possible meron-monopolc 
connection [3j], the monopole solution seems to play a 
special role at the level of the basic structure of the the- 
ory. Multipole solutions have also been discussed in the 
literature j4(. 

Here we want to revisit the monopole solution, address- 
ing the generalized Wu-Yang monopole (WYM) [5] . The 
simple WYM is a Dirac magnetic monopole imbedded in 
the SU(2) gauge group, written in a simple form as: 

A° a = ; A* a = e iab ^ . (1) 

Among recent uses of the WYM, one could cite that the 
description of nuclear motion in a diatomic molecule can 
possess a Wu-Yang like-term [6]. Moreover, the screen- 
ing effects in non-abelian Quark-Gluon Plasmas may give 
rise to Wu-Yang terms as well Q. Furthermore, in Q the 
monopole solution is used in the equation for the inverse 
of the quark propagator. There, a spin-orbit term is com- 
puted together with an energy spectrum. 

Less recent is the belief that a condensed phase of 
monopoles could encode the key of the confinement mech- 
anism. Some authors argue that in the maximal abelian 
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sub-algebra of the Yang-Mills theory, the monopole as- 
sumes a crucial role in long distance effects, such that it 
gives rise to a so-called "Abelian dominance" j^j. 

The WYM has also been used to investigate Gribov 
copies HQ] . This connects with the problem of building a 
suitable generating functional for the Yang-Mills theory 
Green's functions valid beyond perturbation theory. 

In this work, our goal is to obtain solutions of the field 
equations assuming a WYM configuration. With this in 
hands, one can estimate an energy scale and study Gribov 
copies, among other applications. 

The classical configurations that we are going to dis- 
cuss are of the generalized WYM type 



A ia = e iab ^f(r) . (2) 



Of course the generalized WYM keeps its primitive na- 
ture. It mixes the color and the space indices in the 
same way, and includes the original WYM as a particu- 
lar case. The corresponding field equation is a nonlinear 
ordinary differential equation for f(r) that will be inves- 
tigated with a combination of analytical and numerical 
tools. This paper is organized as follows. In section [TT1 
the generalized Wu-Yang configuration is discussed, with 
the computation of F^ v and the classical energy. In IIIII 
the solutions of the classical equation of motion as power 
series is discussed. In section llVl two different numerical 
solutions of the classical field equation are computed and 
their chromomagnetic fields, for large and small r, are 
discussed. In section [V] we explore the conformal invari- 
ance of the theory to introduce a mass scale of the order 
of a typical nonperturbative scale ~ 1 GeV (or ~ 1 fm) 
for the monopole. Then, one can estimate the monopole 
size (energy). Monopole like configurations with discon- 
tinuous first derivatives are discussed in|Vl] Finally, in 
I VIII we resume and conclude. 
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II. GENERALIZED WU-YANG MONOPOLE 
ANSATZ 



In this work, we are going to treat several classical as- 
pects of the Yang-Mills theory, given by the Lagrangian 
density C = — \F^ V F^ U . From now on, it will be as- 
sumed that the rescaling of the gauge fields has been 
performed, so that the coupling constant g will be omit- 
ted. Furthermore, only the SU(2) gauge group will be 
considered, for the sake of simplicity. 

The classical equations of motion are a set of coupled 
non-linear partial differential equations. Among the pos- 
sible classical configurations, one may take the general- 
ized Wu-Yang monopole discussed previously. The time 
component of the gauge field vanishes, Aq = 0, and the 
spatial components are given by equation @ . Note that 
this static configuration satisfies the Lorentz, d^A a ^ = 0, 

and Coulomb , V • A a — 0, gauge fixing conditions auto- 
matically. 

Here one might of course pose the question of the com- 
pleteness of this configuration. Indeed, the f(r) function 
accounts for distance fluctuations. However, in terms of 
directions in the color space, one could think that the 
monopole configuration is very restricted. Further, we 
see that the mixing of color and space-time indices is 
also a quite strong restriction. 

Let us go ahead and calculate the Lagrangian density. 
It is just a matter of algebraic work to find the form of the 
nonabelian gauge tensor due to the generalized Wu-Yang 
configuration. It reads 
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Note that the chromoelectric field vanish. Another im- 
portant object is the energy density, obtained by the con- 
traction of the Fpv tensor. Then, the corresponding chro- 
momagnetic energy density is 
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The total energy is given after integrating §5§ over the 
space volume Now, it is convenient to make the replace- 
ment / = g + 1, which gives the energy as: 
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(6) 



Of course, depending on the form of the integrand, g(r) 
may have a confining shape or not. Note that the require- 
ment of finite energy prevents g(r) of growing without 
bound with r. To achieve the solutions for g(r), one has 
to solve the equations of motion, obtained by requiring 



that g(r) minimizes the energy. The equation of motion 
reads: 



9 



(7) 



This equation, obtained by minimizing the energy, is ex- 
actly the same one that is derived using ([2]) in the classi- 
cal equations of motion of the classical Yang-Mills theory. 
Of course, first one has to solve equation ([7]) to compute 
the classical energy from the generalized WYM configu- 
rations. 



III. SERIES SOLUTION 

Given that ([7]) is nonlinear, it is quite difficult if not 
impossible to get an analytical solution. Then, let us 
construct a series to approximate the solution. First note 
that equation (J7J is symmetric under scale transforma- 
tions: if g(r) is a solution, then g(Xr) is also a solution, 
for any A > 0. To see this, multiply the equation by A 2 
and convert it to an identical equation in the new variable 
x = r/A, whose solution is g(Xx). Next, note that any fi- 
nite solution must remain in the interval — 1 < g(r) < 1. 
If g(r) > 1 with g'(r) > 0, a necessary condition for 
approaching the value 1 from below, the solution grows 
without bound. If g(r) < —1 with g'(r) < 0, a necessary 
condition for approaching —1 from above, the solution 
diminishes without bound. Furthermore, for a solution 
of the equation to be finite in the limits of its domain, 
we must have 



;(r) -+ ± (l - (V) 2 ) 







and 



g(r) -> ± 1 - — 
Afcr 



(8) 



(9) 



for some A; and A^. Further, we note that if g(r) is a 
solution, then —g(r) is as well. 

This is sufficient to give us a hint about the necessary 
form of the required series solution. We cast the solution 
g(r) as: 



g( r ) 



A(r) 



N A 

71 = 



N B 

B{r) = £ b n r- 

n=0 



(0<r<R) 



(R < r) 



(10) 



The reader should be aware that the first terms of these 
series are well known - see, for instance, appendix G of 

Q. 

The function g(r) presents distinct behaviors depend- 
ing on the distance to the origin. For small r, A(r) favors 
growth of the interaction with increasing distance. For r 
greater than R, B(r) favors the decrease of the interac- 
tion with distance. 
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FIG. 1. The series solutions A(r) and B(r). The vertical lines 
are the naive radius of convergence for each series - see text 
for details. 



convergence occurs for r < 1.75, while for the B(r) series 
convergence happens for r > 0.73. 

We have investigated how the naive estimate of the 
radius of convergence changes with the order of the trun- 
cation, i.e. with Na and with Nb- The analysis of the 
series for A(r) and B(r) show that, with the exception 
of the first few terms, the estimate of the radius of con- 
vergence is almost independent of Na and Nb- We can 
then consider matching the series A{r) to the value and 
first derivative of the series B(\r) for some A and some 
r in the interval 0.73/A < r < 1.75. It can be shown, 
however, that such values of A and r do not exist, which 
prevents the construction of a global solution of ([7]) in 
the form of the power series. 

The series solutions can be seen in figure [T] together 
with the naive radius of convergence. 

IV. NUMERICAL SOLUTIONS 



Using the above solution in the differential equation 
(O, a recursion relation for the various a n and another 
for the b n coefficients can be derived. Indeed, to lowest 
order in r or in 1/r for the A(r) and B{r) functions, 
respectively, one has that Og — a = bQ~b = 0. Therefore 
clq = 0, ±1 and bo = 0, ±1. The zero value implies always 
the trivial solution, i.e. a n = and/or b n — 0. Given 
the symmetry g(r) — > —g(r), we will consider only the 
+1 solution. Furthermore, for A{r) it turns out that the 
coefficients for odd powers of r vanish. After redefining 
the A(r) series to include only even powers of r and using 
the +1 value for the clq and bo, the a\ and b\ coefficients 
are arbitrary and the remaining coefficients are given in 
terms of oi and b\. 

In the following, in the computation of Air) and B(r) 
the truncations Na = Nb = 18 are used everywhere. 
For completeness, we write the first terms of the series 
solutions. For a% = b% = 1, it comes that 
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and 
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(12) 



We call the reader attention that, due to conformal in- 
variance, A(r) and A(Xr) or B(r) and B(Xr) are solutions 
of the equations of motion. For the truncation described 
above, one can estimate the radius of convergence of the 
series from the higher order terms in the series. For A(r) 



A numerical solution of equation (J?} can be built com- 
bining the series for A(r) and B{r), described in the 
previous section, with a Taylor expansion of g(r). If at 
r = ro the solution g(r) along with its first and second 
derivatives are known, one can expand g(r) around ro 

g(r + Sr) = g(r ) + g'(r )Sr + lg"(r )(Sr) 2 + ■■■ (13) 

and similar expansions for g'(r) and g"(r) can be writ- 
ten. Replacing, in the differential equation 0, the ex- 
pansion for g(ro + 6r) and g"{ro + <5r), working order by 
order in Sr, one can compute the derivatives g^ n \ro) for 
n = 3, 4, . . . in terms of glr^), g'lr ) and g"{r ). In our 
numerical solutions, discussed below, we have considered 
expansions of glr) up to order 6. 

For the solution called gi n t(f), the integration starts 
at r = 0, setting glr) — A(r), and using an integration 
step of Sr = 10~ 4 . The series solution was used up to 
the point where the relative error between the l.h.s. and 
r.h.s of equation ([7]) was below a certain error, which we 
toke to be 10~ 6 . If the solution g ~ 0, then instead of 
the relative error we have used the absolute error, i.e. the 
absolute difference between the l.h.s. and r.h.s. of the 
differential equation. This procedure (so far, using A(r)) 
gives g(r), g'(r) and g"{r), as the series for A(r) up to 
r = 0.61. From this point onwards, g(r) was propagated 
using the Taylor series procedure described at the begin- 
ning of this section. At each step the relative error, or the 
absolute error, in the differential equation was checked. 
If the precision was above the tolerance demanded, then 
the propagation was restarted using the last point where 
the differential equation was satisfied to calculate the co- 
efficients of the Taylor expansion. The procedure was 
repeated up to the point where \g(r)\ = 1. From this 
point onwards it was assumed that glr) = ±1 for all r, 
with the sign chosen to produce a continuous function. 

For the solution called g ext (r) , integration was started 
at r = 300, setting g(r) = —B(r), and using Sr = — 10~ 4 
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FIG. 2. The solution gint{r) and its derivatives. 




FIG. 3. The solution g ex t(r) and its derivatives. 



as integration step. The procedure to compute g e xt(f) 
follows the same lines as the computation of gi n t (r) , with 
the difference that the solution is now propagated back- 
wards. Deviations from the B(r) series solution happen 
for r < 2.12. 

The numerical integration was checked exploring dif- 
ferent integration steps and the conformal invariance of 
equation (J7J). In particular, the scaling of g'{r) and g"{r) 
under a scale transformation was checked carefully. Fur- 
thermore, relaxing the precision on the differential equa- 
tion, we observed that the numerical solution follow the 
series solution almost up to its naive convergence radius. 
The several tests performed gave us confidence on our 
numerics. 

The curves for <?i ra t(r) and g ex t{ r )i including its first 
and second derivatives, are shown in figures [5] and [31 rc- 




FIG. 4. The energy density for gi n t(r) and g ex t(r). 



spectively. 

Let us discuss first the solution g ex t{r). g e xt(T)\s the 
numerical solution found by Wu and Yang in Al- 
though in figure |3j it seems that the curve touches r = 0, 
this is not the case. A closer look shows that the smallest 
value of r numerically available is 0.0133. Our numeri- 
cal procedure was not able to go below this value of r. 
Anyway g e xt(r) matches perfectly the solution of the lin- 
earized field equation ([7]). i.e. 



g(r) = y/r Go sin 



'V3 



In i 



G\ cos 




(14) 

with Go = -0.079223 and Gi = -0.4251983. g ex t(r) 
oscillates as r approaches zero and its derivative diverges 
at r = 0. The classical energy, as given by (O, associated 
with g e xt{r) is infinite. 

On the other hand, the solution gi n t(r) has finite en- 
ergy, £ = 33.63 in dimensionless units, but its first deriva- 
tive is discontinuous at r = 2.024 where g = — 1. For 
r > 2.024 one assumes a constant solution with g = — 1 
and for r > 2.024 the energy density associated with 
gint(r) vanishes. The energy density for gi n t{r) and 
gextir) is reported in figure 2] 



A. The Chromomagnetic Field 

The chromomagnetic fields are given in terms of / = 
1 + g by equation For each of the numerical solutions 
described previously, one can build two different / = l±g 
and its large distance behavior is given by 



9int{r) 



and 



9ext(r) 



i.e. 



f(r) 




(15) 



(16) 
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and 



9 = +9ext 
9 = ~9ext ■ 



(17) 



For the first solution Fjj = 0. For the second and third 
solutions, F°- ~ 4/r 2 at large r. Finally, for the last solu- 
tion the chromomagnetic field associated with monopole 
vanishes faster and goes like F?- ~ 1/r 3 at large dis- 
tances. Given the behavior of F°j at large distance, one 
can say that solutions two and three are monopole like 
solutions, i.e. F°- ~ 1/r 2 , while the fourth solution is a 
dipole like solution, i.e. Ffi ~ 1/r 3 . 

In what concerns the short distance behavior, for all 
configurations diverges at r = 0. The exception be- 
ing the trivial solution with a vanishing chromomagnetic 
field. 



V. WU-YANG MONOPOLES AND CLASSICAL 
ENERGY 

In the previous section we have computed monopole 
configurations with infinite classical energy (those asso- 
ciated with g e xt{r)), and monopoles with finite energy 
(associated with <?j n t(r)). However, in what concern 
monopoles with finite energy, besides the solutions as- 
sociated with gi n t (r) , the constant solutions of the equa- 
tion of motion ([7]) also belong to this family of functions. 
The equation of motion ([7]) has the constant solutions 
g = 0, ±1 but only 5 = 1 and g = — 1 have finite zero 
energy - see the definition ©. In particular, g = —1 is 
associated with the trivial = configuration and will 
not be considered any further. 

The monopole with g = 0, i.e. / = 1, despite having 
infinite classical energy, has an associated chromomag- 
netic field vanishing as 1/r 2 for large r. 

For the monopole configurations with classical infinite 
energy, the divergence in £ is associated with an ultravi- 
olet, i.e. short distance, divergence of the energy density. 
Like in classical electrodynamics, one can cure the di- 
vergence via the introduction of a short distance cutoff. 
Then, a scale is introduced into the theory and conformal 
invariance is broken explicitly. 

On the other hand, finite energy monopoles are either 
constant solutions of ([7]) or have a discontinuity on the 
first derivative of /. 

Let us go back to the solution gi n t(r). Conformal in- 
variance means that if gi nt (r) is a solution of equation ((7|) , 
5int(Ar) is also a solution of the same differential equa- 
tion with energy 33.63/A. Assuming that this energy is 
a typical hadronic scale ~ 1 GeV, then A = 33.63 GeV -1 
and the discontinuity occurs for r^ = 2.024 A = 13.4 fm, 
roughly one order of magnitude larger than the typical 
hadronic scale ~ 1 fm. If instead, one uses the mass of 
the lowest glueball to set the scale M g i ue =1.7 GeV fl2^ . 
these figures should be divided by 1.7, i.e. A = 19.78 
GeV -1 and r c i — 7.9 fm. Again, the discontinuity on 



the first derivative of / happens at a rather large length 
scale. 

The hadronic length scale ~ 1 fm comes from the anal- 
ysis of baryons and mesons constituted by quarks. Our 
reasoning shows that if the finite energy solution built 
from gint{r) plays a role on the dynamics of the Yang- 
Mills theory, then the large distance behavior is clearly 
dominated by gluon fields. Furthermore, these gluon 
fields have a spatial size about ten times larger than the 
usual quark content. Ignoring, for the moment, possi- 
ble implications of such a picture on quark matter phe- 
nomenology, in the next section we investigate how far 
can one take the discontinuous monopole configurations. 



VI. DISCONTINUOUS WU-YANG 
MONOPOLES WITH FINITE ENERGY 

Discontinuities in the derivative can be introduced in 
the field equation ([7]) via a Dirac-delta function 8(r—R ). 
Away from r = Rq one recovers the original equation and 
g{r) is given by one of the solutions discussed in section 
HI1 Adding a term like ( 8(r — Ro), where £ is a constant, 
to (O is equivalent to add ( g(r)S(r— R ) to the density of 
energy. Remember that the field equation was obtained 
directly from the classical Yang-Mills equations for the 
generalized Wu-Yang configuration and from minimizing 
the energy (|6]). Such a term can be due to the presence 
of static quarks located at Rq. In this sense one can see 
C as the square of the quark wave function at Rq and Rq 
as a mean value of the quark radius. Furthermore, the 
introduction of the new term in £ for a fixed Rq breaks 
conformal invariance. 

A physical scale should be introduced in QCD in order 
to construct the observables. In hadronic physics this is 
set by Aqqd ~ 1 fort -1 , i.e., the scale where conformal 
invariance is broken, as in modern bottom-up approaches 
of the hadronic spectrum within gauge/gravity duality 



(see e.g. [111]). However, our classical solutions, as we 



will see, attain too large energies, in the hadronic scale, 
when the typical extension of the chromagnetic field is 
of the order of 1 fm, which suggests that it is required a 
collective phenomena to excite such field configurations. 
More on that discussion will come later on. 

So let us assume that the classical field equation is 
given by 
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= 4( 5 2 -l) + (5(r-R 



(18) 



For r ^ R the solutions of (fTB")) are the functions g(r) 
discussed in section [II]- see figure [S] Assuming that g{r) 
is continuous at Rq, integrating (IT8t from Rq — e to Ro + e 
gives, in the limit where e — > + , 



g'(R + e)-g'(R -e) = t. 



(19) 



Finite energy solutions require g(0) = ±1. Therefore, 
for r < Rq only ±gi nt (r) and the constant solutions g = 
±1 are acceptable. On the other side, for r > Rq, there 
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FIG. 5. The solutions gint{r) and g £X t(r) 







■ 




— tfi) 






- gjW 

- gjW 
■-■ g,M 










\ 

\ 

\ 

\ 

\ 

v_ 


\ 

\ 






\ 






\ 




, , , 



FIG. 6. The solutions gi(r) to g±(r). 



is no reason to exclude any of the solutions ±<?i„ t (r), 
igext(f)i g( r ) = or g(r) = ±1 which can be combine 
to produce several finite energy solutions. The solution of 
equation (fT5|) gives an / continuous and /' discontinuous 
at Rq - see equation (fT9|) . The various solutions are 



(20) 




g ln t(r) r<R Q = 1.1831 







r>R 



93 (r) = 



g in t{r) r<R = 2.0240 

-1 r > R 

g in t{r) r < i? - 0.8651 
-g e xt(r) r > R 



(21) 



(22) 



with a (3 = 1.5432 and 

g m t(r) r < R = 1.6685 
54 M = < (23) 

with (2 = 1-009. In all cases C was computed with equa- 
tion (fTT)|) . The various solutions can be seen in figure 

m 

Assuming that the energy density is modified by the 
term £ g(r) S(r — Rq), then the total energy is given by 



«+oo 



£ = 4tt 



dr (g'f + 



-1 



2r 2 



-1, 53(^0) 



(g(Ro)J ■ 

(24) 
0.3859 and 



Since gi(R ) = 0, 52(^0) 
54(i?o) = —0.5877, the energy of the monopoles is, in 
dimensionless units, £ 1 = 17.28, £ 2 = 19.18, £ 3 = 20.90, 
£4 = 21.06, respectively. In physical units r — Xf, £ = 
£/X, where the quantities with tilde are dimensionless. 
Setting the energy of the monopole to 1 GeV, then Ai = 
17.28 GeV" 1 = 3.410 fm, A 2 = 19.18 GeV" 1 = 3.785 
fm, A 3 = 20.90 GeV- 1 = 4.124 fm, A 4 = 21.06 GeV" 1 = 
4.156 fm. 

In physical units, the discontinuity in the derivative 
happens at XRq. For gi(r), this means 4.0 fm. For g2(r) 
one gets A i?o = 7.7 fm. For gs(r) one gets XRq = 3.6 fm 
and for g±{r) it follows that XRq = 6.9 fm. Note that if 
one uses the lowest glueball mass 1.7 GeV [l2| to set the 
scale, one should divide these figures by 1.7, i.e. the dis- 
continuities in the derivative would occur at 2.4 fm, 4.5 
fm, 2.1 fm, 4.1 fm for gi(r), 9z(r), gz{r). g 4 {r), respec- 
tively, beyond which the chromomagnetic fields would de- 
cay, at least, as 1/r 2 . Again, the scale at which = Ro 
happens is larger, at least two times, than the typical 
hadronic scale. 

In dimensionless units, the energy density for the var- 
ious solutions under discussion is reported in figure [7] 
Another measure of the monopole size is given by the 
maximum in the energy density. For gi[r), gi (r) and 
g±(r) the maximum occur at r = 1. For g% (r) the maxi- 
mum occur at r = 0.8651. Given the A values computed 
in the previous paragraph, if one assumes a scale of 1 
GeV, then the monopole sizes are 3.4 fm, 4.5 fm, 3.6fm 
and 4.2 fm for gi(r), gzfr), gz{r) and g±{r), respectively. 
Otherwise if one uses the glueball mass to set the scale, 
the characteristic sizes are 2.0 fm, 2.2 fm, 2.1 fm and 2.4 
fm, respectively. 



VII. RESULTS AND DISCUSSION 

In this work we have revisited the Wu-Yang monopole 
for pure SU(2) Yang-Mills theory. Assuming that the 
function f(r) = 1 + g(r) in equation ([2]) can be written 
as a power series defines, together with the requirement 
of finite energy, the boundary conditions at infinity and 
at r — 0, namely: g(Q) — ±1 and 5(00) = ±1. Besides 
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FIG. 7. Energy density, in dimensionless units, for the fi- 
nite energy solutions gi(r), gi(r) and gz{r). Note the 1/r 2 
dependence at large r for gi(r), 53 (r) and g± (r). 



the constant solutions g(r) = and g{r) — 1, Wu and 
Yang were able to find, long ago, a numerical solution 
with (7(0) = and g{oo) = —1. 

In order to solve the classical equations of motion, we 
rederived the asymptotic series solutions, close to r = 0, 
named A(r), and for large r, named B(r). The first terms 
of the series solutions are given in equations (jlip and 
(fT2"f . respectively. We have investigated the two series 
and tried to build a solution of the equation of motion 
matching A{r) up to order r 18 and B(r) up to order 
1/r 18 . It turns out that the radius of convergence make 
such a type of solution impossible. Indeed, we found 
that increasing the order of the expansions leaves almost 
unchanged the naive radius of convergence discussed at 
the end of section IIII1 

In section IIVI numerical solutions of the equations of 
motion were computed and the behavior of their chro- 
momagnetic fields discussed at small and large distances. 
Our numerical procedure found two different solutions. 
The numerical solution with g(0) = and g(oo) = — 1 
discovered by Wu-Yang, which we named gextir), and a 
new solution <?mt(r). 

In what concerns g e xt{f)i we were able to provide an 
analytical description of its behavior at small distances 
in (fH)) . Further, the classical energy associated with 
gext{r) diverges due to its behavior at r ~ 0. On the 
other hand, the energy associated with gi n t(r) is finite 
but one is forced to introduce a discontinuity on the first 
derivative of f(r). Anyway, we have explored the con- 
formal invariance of the theory to set a scale for gi nt (r) 
assuming that its classical energy is a typical hadronic 
scale. It turns out that the discontinuity, assumed to 
give a measure of the monopole size, should happen at a 
length scale which is, at least, four times larger than the 
typical hadron size. Of course, if one requires that the 
discontinuity reproduces the typical hadronic size, then 
is the monopole energy which is to large, at least by a 
factor of 8, than a typical hadronic mass. However, if 



one uses the maximum of the energy density to define 
the monopole spatial extent, we get a size of ~ 6.6 fm 
or ~ 3.9 fm assuming a classical energy of 1 GeV or 1.7 
GeV, respectively, for the monopole. 

In section [VI] we explore solutions of the classical equa- 
tion of motion which allow for discontinuities on the first 
derivative. In particular, a Dirac delta function is added 
to equation (J7J. As discussed, such a contribution can 
be motivated by the presence of quarks. The introduc- 
tion of a Dirac delta function 5(r — Rq) with a fixed i?o, 
breaks explicitly invariance of the equations of motion 
under a scale transformation. Three different discontinu- 
ous monopoles configurations (j2"0")) - (|2"2"|) were computed. 
If one uses a typical hadronic scale to define the monopole 
energy, then, depending on the scale, Rq = 7.7 fm or 
i?o = 2.4 fm or larger. Again, this is quite a large length 
scale. Probably, this means that the monopole configu- 
rations computed here play a minor role, if any, on low 
energy phenomenology. 

Curiously, if one takes the critical energy density e c = 
0.7 GeV/fm 3 for the transition to the plasma phase [l3[ 
to set the scale, i.e. one puts e c — (£/A)/(4-7ri? 3 /3), where 
E is the dimensionless monopole energy, and take R = 1 
fm, a typical hadronic length, then one gets a monopole 
energy of 2.9 GeV. Using again the discontinuity on the 
first derivative as a measure of the monopole size, it fol- 
lows that for the solution <7» n t(r), described in section 
IIVI A = 2.26 fm and Ro — 4.6 fm. Remember that in 
section IIVI Ro was estimated to be 13.4 fm or 7.9 fm, 
depending on the choice of scale. For the solution <?i(r) 
of section fVTl A = 1.16 fm and Rq — 1.38 fm, to be com- 
pared with 4.0 fm or 2.4 fm found in section I VII For 
the solution g2{r), A = 1.29 fm and Rq = 2.6 fm, to be 
compared with 7.7 fm or 4.5 fm found in section IVT1 For 
the solution <73(r), A = 1.41 fm and i?o = 1-2 fm, to be 
compared with 3.6 fm or 2.1 fm found in section IVT1 For 
the solution <74(r), A = 1.42 fm and i?o = 2.4 fm, to be 
compared with 6.9 fm or 4.1 fm found in section lVTl This 
suggests that these configurations may play an important 
role in the confining region and in the transition to the 
plasma phase. Furthermore, the so-called discontinuous 
solutions, with sizes in the range 1.2 — 2.6 fm, are con- 
siderable shorter than gi n t{r). The difference is certainly 
due to the Dirac-delta term introduced in the equation 
and can be seen as a clear sign of the role of the break- 
ing of conformal invariance due to the coupling to the 
quark fields. Hopefully, one would be able to investigate 
possible signatures for identifying the monopole configu- 
rations in the transition to the plasma phase, in ongoing 
and future experiments to probe this state of matter with 
high-energy heavy ion collisions. 
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